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This paper discusses the IAU Resolutions B1.3, B1.4, B1.5 and B1.9 (2000) that were 
adopted during the 24th General Assembly in Manchester, 2000 and provides details 
and explanations for these Resolutions. It is explained why they present significant 
' progress over the corresponding IAU 1991 Resolutions and why they are necessary in 

the light of present accuracies in astrometry, celestial mechanics and metrology. In fact 
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most of these Resolutions are consistent with astronomical models and software already 
in use. 

The metric tensors and gravitational potentials of both the Barycentric Celestial 
Reference System and Geocentric Celestial Reference System are defined and discussed. 
The necessity and relevance of the two celestial reference systems are explained. The 
transformations of coordinates and gravitational potentials are discussed. Potential 
coefficients parameterizing the post-Newtonian gravitational potentials are expounded. 
Simplified versions of the time transformations suitable for modern clock accuracies are 
elucidated. Various approximations used in the Resolutions are explicated and justified. 
Some models (e.g. for higher spin moments) that serve the purpose for estimating orders 
of magnitude have actually never been published before. 

Subject headings: relativity, astrometry, celestial mechanics, reference systems, time 

1. Introduction 

It is clear that, beyond some threshold of accuracy, any astronomical problem has to be 
formulated within the framework of Einstein's theory of gravity (General Relativity Theory, GRT). 
Many high precision astronomical techniques have already passed this threshold. For example, 
Lunar Laser Ranging (LLR) measures the distance to the Moon with a precision of a few cm, thereby 
operating at the 1CT 10 level. At this level several relativity effects are significant and observable. 
Relativity effects related with the motion of the Earth-Moon system about the Sun are of the order 
(forbitai/ c ) 2 — 10 -8 . The Lorentz contraction of the lunar orbit about the Earth that appears in 
barycentric coordinates has an amplitude of about 100 cm whereas in some suitably chosen (local) 
coordinate system that moves with the Earth-Moon barycenter the dominant relativistic range 
oscillation reduces to a few cm only (Mashhoon 1985; Soffel et al. 1986). 

The situation is even more critical in the field of astrometry. It is well known that the grav- 
itational light deflection at the limb of the Sun amounts to 1.75" and decreases only as 1/r with 
increasing impact parameter r of a light ray to the solar center. Thus, for light rays incident at 
about 90° from the Sun the angle of light deflection still amounts to 4 mas. To describe the accu- 
racy of astrometric measurements it is useful to make use of the parameter 7 of the parameterized 
post-Newtonian (PPN) formalism. We would like to emphasize that this paper deals solely with 
Einstein's theory of gravity, where 7 = 1, and not with the PPN formalism. Nevertheless, the in- 
troduction of 7 is useful if one talks about measuring accuracies. In the PPN formalism, the angle 
of light deflection is proportional to (7 + l)/2 so that astrometric measurements might be used for 
a precise determination of the parameter 7. Meanwhile, VLBI has achieved accuracies of better 
than 0.1 mas, and regular geodetic VLBI measurements have frequently been used to determine 
the space curvature parameter. A detailed analysis of VLBI data from the projects MERIT and 
POLARIS/IRIS gave 7 = 1.000 ± 0.003, (Robertson et al. 1984; Carter et al. 1985) where a 
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formal standard error is given. Recently an advanced processing of VLBI data provided the best 
current estimates 7 = 0.9996 ± 0.0017 (Lebach et al. 1995) and 7 = 0.99994 ± 0.00031 (Eubanks 
et al. 1997). Current accuracy of modern optical astrometry as represented by the HIPPARCOS 
catalog is about 1 milliarcsecond, which gave a determination of 7 at the level of 0.997 ± 0.003 
(Froeschle, Mignard & Arenou 1997). Future astrometric missions such as SIM and especially 
G AIA will push the accuracy to the level of a few microarcseconds (//as) , and the expected accu- 
racy of determinations of 7 will be 10 -6 — 10 -7 . The accuracy of 1 /tas should be compared with the 
maximal possible light deflection due to various parts of the gravitation field: the post-Newtonian 
effect of 1.75" due to the mass of the Sun, 240 //as caused by the oblateness of Jupiter, J2 (10 //as 
due to Jupiter's J4), the post-post-Newtonian effect of 11 //as due to the Sun, etc. This illustrates 
how complicated the relativistic modeling of future astrometric observations will be. It is clear that 
for such high accuracy the corresponding model must be formulated in a self-consistent relativistic 
framework. 

Another problem worth mentioning is that of time measurement. The realization of the SI 
second (the unit of proper time) has improved by one order of magnitude in the last few years with 
the advent of laser cooled atomic clocks (Lemonde et al. 2001; Weyers et al. 2001, and references 
therein), and is now below two parts in 10 15 . This should be compared with the dimensionless 
quantity Ue/c 2 ~ 7 x 10~ 10 , which gives the order of magnitude of relativistic effects produced by 
the gravity field of the Earth itself in the vicinity of its surface. In the near future, laser cooled 
atomic clocks in micro-gravity are expected to lead to a further improvement by at least one order 
of magnitude. At present several clock experiments in terrestrial orbit are planned, such as the 
Atomic Clock Ensemble in Space project (ACES, Lemonde et al. 2001). These in turn are likely 
to lead to clock experiments in solar orbits like the Solar Orbit Relativity Test (SORT) project. 
All of these experiments require a detailed account of many subtle relativistic effects. 

Finally we would like to mention the problem of geodesic precession-nutation (a relativistic 
effect that is discussed in more detail below, Misner, Thorne & Wheeler 1973; Soffel 1989) and the 
description of Earth's rotation in a suitably chosen Geocentric Celestial Reference System (GCRS). 
Geodesic precession amounts to 1.9" /century and geodesic nutation is dominated by an annual term 
with amplitude 0.15 mas. Since the GCRS is chosen to be kinematically non-rotating geodesic 
precession-nutation should be contained in the model describing the relation between the GCRS 
and the International Terrestrial Reference System (ITRS). According to the IAU 2000 Resolution 
B1.6 this relativistic precession- nutation is indeed contained in the present IAU precession-nutation 
model. 

These examples show clearly that high-precision modern astronomical observations can no 
longer be described by Newtonian theory, but require Einstein's theory of gravity. 

The consequences of this are profound for the basic formalism to be used since one often tends 
to express it in terms of "small relativistic corrections" to Newtonian theory. This can lead to 
misconceptions and mistakes. One central point is that in Newton's theory, globally preferred co- 
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ordinate systems exist that have a direct physical meaning. In the Newtonian framework, idealized 
clocks show absolute time everywhere in the universe at all times, and global spatial inertial co- 
ordinates exist in which dynamical equations of motion show no inertial forces. This is no longer 
true in GRT. Usually space-time coordinates have no direct physical meaning and it is essential 
to construct the observables as coordinate independent quantities, i.e., scalars, in mathematical 
language. This construction usually occurs in two steps: first one formulates a coordinate picture 
of the measurement procedure and then one derives the observable out of it. This leads us to 
the problem of defining useful and adequate coordinate systems in astronomy. The underlying 
concept in relativistic modeling of astronomical observations is a relativistic four-dimensional ref- 
erence system (RS). By reference system we mean a purely mathematical construction (a chart or 
a coordinate system) giving "names" to space-time events. In contrast to this a reference frame 
is some materialization of a reference system. In astronomy the materialization is usually realized 
by means of a catalog or ephemeris, containing positions of some celestial objects relative to the 
reference system under consideration. Hence it is very important to understand that any reference 
frame is defined only through a well-defined reference system, which has been used to construct 
physical models of observations. 

In the following, a 4-dimensional space-time reference system will be described by four coor- 
dinates x a = (x°,x l ) = (x° . X • X - X 3 ). Here and below the greek indices (e.g., a) take the values 
0, 1, 2, 3 and the latin ones (e.g., i) take the values 1, 2, 3. The index refers to the time variable 
and the indices 1,2,3 refer to the three spatial coordinates. For dimensional reasons one usually 
writes x° = ct where c is the speed of light and t a time variable. According to the mathematical 
formalism of General Relativity, a particular reference system is fixed by the specific form of the 
metric tensor g a/ 3(t,x l ), which allows one to compute the 4-distance ds between any two events x a 
and x a + dx a according to the rule 

ds 2 = g a p{t, x 1 ) dx a dx 13 = gooc 2 dt 2 + Ig^c dt dx % + gijdx l dx^ , (1) 

where Einstein's summation convention (summation over two equal indices) is implied. The metric 
tensor allows one to derive translational and rotational equations of motion of bodies, to describe 
the propagation of light, and to model the process of observation. Examples of such modeling 
include relating the observed (proper) time of an observer to the coordinate time t, or relating the 
angles between two incident light rays as observed by that observer to the corresponding coordinate 
directions. All of these components can be combined into a single relativistic model for a particular 
kind of observations. Such a model contains a certain set of parameters describing various properties 
of the objects participating in the process of observations. These parameters should be determined 
from observations. Many of these parameters crucially depend on the reference system used to 
formulate the model of observations (e.g., the initial positions and velocities of certain bodies). 
Some other parameters might not depend at all on the reference system (e.g., the speed of light in 
vacuum). On the other hand, according to the principle of covariance different reference systems 
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covering the region of space-time under consideration are mathematically equivalent in the sense 
that any such system can be used to model the observations. This freedom to choose the reference 
system can be used to simplify the models or to make the resulting parameters more physically 
adequate. 

It is widely accepted that in order to describe adequately modern astronomical observations one 
has to use several relativistic reference systems. The solar system Barycentric Celestial Reference 
System (BCRS) can be used to model light propagation from distant celestial objects as well as the 
motion of bodies within the solar system. The Geocentric Celestial Reference System is physically 
adequate to describe processes occurring in the vicinity of the Earth (Earth's rotation, motion 
of Earth's satellites etc.). The introduction of further local systems (selenocentric, martian etc.) 
might be of relevance for specific applications, where physical phenomena in the vecinity of the 
corresponding body play a role. 

The necessity to use several reference systems can be understood from the following example. 
If we were to characterize terrestrial observers by the difference between their BCRS coordinates 
and the BCRS coordinates of the geocenter, the positions of the observers relative to the geocenter 
would be altered by time-dependent, relativistic coordinate effects (such as Lorentz contraction) 
which have nothing to do with the Earth's rotation or with geophysical factors and which would 
vanish if one employs suitable GCRS coordinates instead. On the other hand, the coordinate 
positions derived with VLBI observations are used to investigate local geophysical processes and 
some adequate geocentric reference system allows one to simplify their description. 

The basic idea is to construct a special local reference system for each material subsystem, 
where relativistic equations of motion of a test body inside the considered subsystem take a partic- 
ularly simple form. In such a local reference system the influence of external matter, in accordance 
with the equivalence principle, should be given by tidal potentials only, that is, by potentials whose 
expansions in powers of local spatial coordinates in the vicinity of the origin of the corresponding 
reference system starts with the second order (the linear terms representing inertial forces may also 
exist, but can be eliminated if desired by a suitable choice of the origin of the local coordinates). 

Two advanced relativistic formalisms have been worked out to tackle this problem in the 
first post-Newtonian approximation of General Relativity. One formalism is due to Brumberg and 
Kopeikin (Brumberg & Kopeikin 1989; Kopeikin 1988, 1990; Brumberg 1991; Klioner & Voinov 

1993) and another one is due to Damour, Soffel and Xu (Damour, Soffel & Xu 1991, 1992, 1993, 

1994) , frequently called the DSX framework. The IAU 2000 Resolutions B1.3-B1.5 are based on 
these approaches. These Resolutions extend corresponding older ones that are reconsidered in the 
next section. From a mathematical point of view Resolution B1.3 recommends the use of certain 
coordinates and the way of writing the metric tensor. Clearly one might use any coordinate system 
that might be well adapted to the specific problem of interest. Nevertheless because of the high 
risk of possible confusion the strategy of recommending special coordinate systems (to fix the gauge 
completely in the mathematical language) has significant advantages. If different coordinates arc 
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employed to derive certain results this should be stressed explicitly so that they can be transformed 
into the reference systems recommended by the IAU and can be compared with the results derived 
in the IAU framework. 

The organisation of the present paper is as follows (see also Figure 1). In Section 2 the prin- 
cipal content of the IAU 1991 recommendations on relativity (Section 2.1) and the further related 
IAU and IUGG resolutions (Section 2.2) are repeated. The IAU 2000 resolutions on relativity 
(Resolutions B1.3, B1.4, B1.5 and B1.9) are discussed in Section 3. The full text of the IAU 2000 
resolutions on relativity is given in Appendix A. Section 3.1 briefly clarifies the necessity and the 
role of the two celestial reference systems defined by the IAU resolutions. The Barycentric Celestial 
Reference System (BCRS) defined by Resolution B1.3 is discussed in Section 3.2. Section 3.3 is de- 
voted to a discussion of the Geocentric Celestial Reference System (GCRS) as well as the definition 
of the geocentric gravitational potentials also defined by Resolution B1.3. The coordinate trans- 
formations between the BCRS and GCRS, also fixed by Resolution B1.3, are explained in Section 
3.4. Potential coefficients which can be used to represent in a meaningful way the post-Newtonian 
geocentric gravitational potential of the Earth in its immediate vicinity are fixed by Resolution 
B1.4 and explained in Section 3.5. As an illustration the gravitational potentials of the BCRS are 
calculated in Section 3.6 for the simplified case when all gravitating bodies of the solar system can 
be characterized by their masses only (no further structure of the gravitational field of the bodies 
is considered). A similar form of the barycentric gravitational potentials is used in Resolution B1.5 
where a practical relativistic framework for time and frequency applications in the solar system is 
formulated. This practical relativistic framework is discussed in Section 4. The practical transfor- 
mation between the coordinate times of the BCRS and GCRS is explained in Section 4.1 while the 
transformations between various kinds of the time scales appropriate for the Earth's vicinity are 
discussed in Section 4.2. Appendix B contains an explicit proof that the BCRS metric coincides 
with well-known results from the literature. 



The IAU Resolution A4 (1991) contains nine recommendations, the first five of which are 
directly relevant to our discussion. 

In the first recommendation, the metric tensor for space-time coordinate systems (t, x) centered 
at the barycenter of an ensemble of masses is recommended in the form 



2. The IAU 1991 framework and previous recommendations 



2.1. 



The IAU 1991 recommendations 



5oo = -1 + 



2U(t,x) 



Barycentric Celestial RS: 

Coordinate time: t = TCB 

Spatial coordinates: x' 

Metric tensor: g IJV 

Gravitational potentials: w, w' 



Resolution B1.3 
Section 3.2 



Coordinate Transformations: 

fax') < — y (T,X a ) 



Resolution B1.3 
Section 3.4 



Geocentric Celestial RS: 

Coordinate time: T = TCG 

Spatial coordinates: X" 

Metric tensor: G a s 

Gravitational potentials: W, W" 



Resolution B1.3 
Section 3.3 



Transformation of the Potentials: 

(ujX) < — y (W,W a ) 



Resolution B1.3 
Section 3.3 



Geocentric multipole moments: 

M Ll S L 



Resolution B1.4 
Section 3.5 



Practical Transformation: 

TCB < — > TCG 



Resolution B1.5 
Section 4.1 



Fig. 1. — Notations used for quantities of the Barycentric and Geocentric Celestial Reference 
Systems (coordinates, metric, potentials, and multipole moments) with references to the Sections 
and Resolutions where they appear. 
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g 0l = 0(c~ 3 ), (2) 

9ij=Sij (l + ^x^) + 0{c-*) , 

where c is the speed of light in vacuum (c = 299792458 m/s) and U is the Newtonian gravitational 
potential (here a sum of the gravitational potentials of the ensemble of masses, and of an external 
potential generated by bodies external to the ensemble, the latter potential vanishing at the origin). 
The algebraic sign of U is taken to be positive and it satisfies the Poisson equation 



AU = -4ttG P . (3) 

Here, G is the gravitational constant, p is the matter density and A is the usual Laplacian A = 
d 2 /dx 2 + d 2 /dy 2 + d 2 /dz 2 , x = (x,y, z). This recommendation recognizes also that space-time 
cannot be described by a single coordinate system. The recommended form of the metric tensor 
can be used not only to describe the barycentric reference system of the whole solar system, but also 
to define the geocentric reference system centered in the center of mass of the Earth with a suitable 
function instead of U, now depending upon geocentric coordinates. In analogy to the geocentric 
reference system a corresponding reference system can be constructed for any other body of the 
solar system. 

In the second recommendation, the origin and orientation of the spatial coordinate grids for 
the barycentric and geocentric reference systems are defined. Notably it is specified that the 
spatial coordinates of these systems should show no global rotation with respect to a set of distant 
extragalactic objects. It also specifies that the SI (International System of Units) second and the 
SI meter should be the units of time and length in all coordinate systems. It states in addition that 
the time coordinates should be derived from an Earth atomic time scale. 

The third recommendation defines TCB (Barycentric Coordinate Time) and TCG (Geocentric 
Coordinate Time) as the time coordinates of the BCRS and GCRS, respectively. Here we write 
(t = TCB,x % ) and (T = TCG,X % ) for the respective coordinates. The recommendation also 
defines the origin of the times scales in terms of International Atomic Time (TAI) . The reading of 
the coordinate time scales on 1977 January 1, h m s TAI {JD = 2443144.5 TAI) must be 1977 
January 1, h m 32.184 s ). Finally, the recommendation declares that the units of measurements 
of the coordinate times of all reference systems are chosen so that they are consistent with the 
SI second. The relationship between TCB and TCG is then given by the time part of the full 
4-dimensional transformation between the barycentric and geocentric reference systems 



TCB - TCG = c 



t /„,2 



+ C/ext(x E ) dt + ^ E r% 



+ 0{c-% 



(4) 



where x l E and v E are the barycentric coordinate position and velocity of the geocenter, r l E = x l — x\ 
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with x l the barycentric position of some observer, and f7 e xt(x£) is the Newtonian potential of all 
solar system bodies apart from the Earth evaluated at the geocenter. 

In the fourth recommendation another coordinate time, Terrestrial Time (TT), is defined. It 
differs from TCG by a constant rate only 

TCG — TT = Lq x (JD — 2443144.5) x 86400, Lq ~ 6.969291 x 10~ 10 , (5) 

where JD is TAI measured in Julian days, so that the mean rate of TT agrees with that of 
the proper time of an observer situated on the geoid up to a certain accuracy limit. Up to a 
constant shift of 32.184 s, TT represents an ideal form of TAI, the divergence between them being 
a consequence of the physical defects of atomic clocks. It is also recognized that TT is nothing but 
a rescaling of the geocentric coordinate time TCG. 

The fifth recommendation states that the old barycentric time TDB may still be used where 
discontinuity with previous work is deemed to be undesirable. Let us note, however, that TDB was 
never defined in a self-consistent and exact manner. For that reason it cannot be used in theoretical 
considerations. In the notes to the third recommendation the relation of the TCB with TDB is 
given as 

TCB -TDB = L B x {JD - 2443144.5) x 86400, Lb ~ 1.550505 x 10~ 8 . (6) 

Note, however, according to the IAU Resolution C7 (see Section 2.2) JD is defined in Terrestrial 
Time (TT) which makes this formula problematic. 

2.2. Further resolutions 

Resolution 2 (1991) of the International Union of Geodesy and Geophysics (IUGG) defined a 
Conventional Terrestrial Reference System (CTRS) as a reference system resulting from a [time- 
dependent] spatial rotation of the geocentric reference system defined by the IAU 1991 Recommen- 
dations, the spatial rotation being chosen such that the CTRS has no global residual rotation with 
respect to horizontal motions at the Earth's surface. The coordinate time of the CTRS coincides 
with TCG. 

IAU Resolution C7 (1994) recommends that the epoch J2000 as well as the Julian (ephemeris) 
day are defined in TT. 

IAU Resolution B6 (1997) has supplemented the framework by one more recommendation 
stating that no scaling of spatial axes should be applied in any reference system (even if a scaled 
time coordinate like TT is used). Note, however, that this Resolution has been ignored in the 
construction of the International Terrestrial Reference Frame (ITRF) which is defined not with the 
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GCRS spatial coordinates X but with scaled coordinates Xtt = (1 — Lq) X. 

3. The IAU 2000 Resolutions on relativity 

The IAU 1991 framework is unsatisfactory from many points of view. The Einstein-Infeld- 
Hoffmann equations of motion which have been used since the 70s to construct the JPL numerical 
ephemerides of planetary motion cannot be derived from the metric (2). In other words for the 
motion of massive solar system bodies metric (2) is not the post-Newtonian metric of Einstein's 
theory of gravity In the years prior to the XXIII General Assembly in Kyoto (1997) it became 
obvious that the IAU 1991 set of recommendations concerning relativity in astrometry, celestial 
mechanics and metrology is not sufficient for achievable accuracies. Especially with respect to 
planned astrometric missions with microarcsecond accuracies extended and improved resolutions 
became indispensable. For that reason the IAU Working Group "Relativity for astrometry and 
celestial mechanics" together with the BIPM/IAU Joint Committee on relativity for space-time 
reference systems and metrology suggested such an extended set of Resolutions (B1.3 - B1.5 and 
B1.9) that was finally adopted at the IAU General Assembly in Manchester in the year 2000. The 
relevant Resolutions can be found in Appendix A. It is clear that because of their brevity they need 
additional explanations and one has to show how they work in practice. This paper now presents 
a detailed explanatory supplement for these IAU 2000 Resolutions. 

3.1. The role of the two Celestial Reference Systems, BCRS and GCRS 

Some of the reasons why two different celestial astronomical reference systems have to be 
introduced have already been mentioned in the introduction. Here we would like to deepen this 
discussion in several respects. It is clear that for many applications in the fields of astrometry, 
celestial mechanics, geodynamics, geodesy etc. some quasi-inertial or 'space-fixed' reference system 
has to be introduced. Resolution B1.3 defines actually two different Celestial Reference Systems: 
the Barycentric Celestial Reference System (BCRS) and the Geocentric Celestial Reference System 
(GCRS). 

In Newtonian theory one can easily introduce inertial space-time coordinates that cover the 
entire universe. Such inertial coordinates in Newton's theory are unique up to the choice of origin, 
scales, the orientation of spatial axes and up to a constant velocity of origin. In astronomy concep- 
tually we may talk about two different relevant celestial systems: a barycentric and a geocentric 
one that basically serve different purposes. The barycentric celestial system is considered to be 
inertial (external galactic and extragalactic matter being normally neglected) and is used for solar 
system ephemerides, for concepts such as an ecliptic, for interplanetary spacecraft navigation etc. 
The positions of remote objects can be defined in that system. The barycentric celestial system 
presents the fundamental astrometric system where concepts as 'proper motion' or 'radial velocity' 
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can be defined. 

On the other side the geocentric celestial system might be called quasi-inertial since the spatial 
axes are non-rotating in the Newtonian absolute sense whereas the geocenter is accelerated. It is 
employed for the description of physical processes in the vicinity of the Earth, for satellite theory, 
the dynamics of the Earth including Earth's rotation, etc. It is also used for the introduction of 
concepts like equator and the ITRS. Let us denote the time and space coordinates of the barycentric 
celestial system by (t, x), those of the geocentric celestial system by (T, X). In Newton's framework 
the relation between these two sets of coordinates is trivial 

T = t, X = x-x B (t), 

where xg(t) denotes the barycentric position of the geocenter. Because these relations are so 
trivial for some purposes the barycentric and the geocentric celestial systems are not always clearly 
distinguished in the Newtonian framework. 

Of course for astrometric problems one always distinguished between the two celestial systems 
and apparent places of stars from true (barycentric) places. However, annual parallax and aberra- 
tion were merely understood as correction terms that have to be applied to get the 'true' positions 
for the realization of the astronomical quasi-inertial, space-fixed celestial system. Note that the 
definition of the classical astronomical (a, S) system uses concepts from both systems: some eclip- 
tic from the barycentric celestial system and some Earth's rotation pole, Celestial Ephemeris Pole 
(CEP) or Celestial Intermediate Pole (CIP), and its corresponding equator from the geocentric 
celestial system. 

In Relativity theory the situation is actually more complicated. Even in the absence of grav- 
itational fields and a uniformly moving geocenter the two coordinate systems are related by a 
4-dimensional Lorentz transformation of Special Relativity. In our solar system BCRS and GCRS 
coordinates are related by a complicated 4-dimensional space-time transformation (a generalized 
Lorentz transformation) that also contains acceleration terms and gravitational potentials. This 
implies that the two astronomical Reference Systems, the BCRS and the GCRS are actually quite 
different. This has profound consequences for a lot of classical astronomical concepts. 

The BCRS is the basic astrometric celestial reference system. Usually one considers the solar 
system to be isolated, i.e., one ignores all matter and fields outside the system and assumes the 
gravitational potentials to vanish far from the system. It is obvious that ignoring the Galaxy and 
extragalactic objects is an unphysical idealization for several specific questions that, however, will 
not be touched here. If the solar system is considered to be isolated we might follow light-rays 
from some very remote source back in time to the region |x| — > oo that might be called the celestial 
sphere. In the vicinity of the celestial sphere a certain light ray defines spherical angles that might 
appear as catalog values. Actually for reference stars the physical distance from the Earth usually 
plays a role. In that case we might associate any star with a corresponding BCRS coordinate 
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position that will be a function of TCB. From this position vector spherical angles (a*,<5*) can 
be introduced in a very simple manner by 



that can be considered as catalog values. If the coordinate distance of some source tends to infinity 
the two constructions for an astrometric position will coincide. From x*(i) quantities such as 
'proper motion' or 'radial velocity' can be defined as coordinate quantities in the BCRS. Note that 
the problem of 'radial velocity' has exhaustively been discussed by Lindegren & Dravins (2003) 
(see also IAU 2000 Resolutions CI and C2 (Rickman 2001)). Other fields of application of the 
BCRS are solar system ephemerides, interplanetary navigation, etc. 

The definitions of the BCRS given by the IAU 2000 Resolution B1.3 do not fix the orientation 
of spatial axes uniquely but only up to some constant, time independent rotation matrix about the 
origin. One natural choice of orientation is provided by the International Celestial Reference System 
(ICRS). Actually for the construction of the International Celestial Reference Frame (ICRF) and its 
optical counterpart, the Hipparcos catalog, the recommended form of the barycentric metric tensor 
has already been used explicitly in the underlying models. This implies that a set of definitions 
that fix the ICRS completely contains the BCRS definitions. 

There might be other useful possibilities for the orientation of barycentric spatial coordinates. 
One possibility is the orientation according to some ecliptic £q at a certain epoch to defined by 
corresponding solar system ephemerides. Such an ecliptic would coincide with the x-y plane of a 
BCRS [.So] that might be useful for reasons of historical continuity. 

On the other hand quantities and concepts related with the physics in the immediate vicinity 
of the Earth should be formulated in the GCRS. This concerns the gravity field of the Earth 
itself, satellite theory and especially applies to theories of Earth's rotation and their parameters. 
Clearly the spatial GCRS coordinates X can be used to define corresponding unit vectors at the 
geocenter that might be employed to compute spherical angles («gcrs > ^GCRS ) which might be 
called 'geocentric places'. Note, however, that the coordinates of the remote astronomical sources 
are defined in the BCRS only. The calculated GCRS places («GCRS) ^gcrs) are determined by 
incident light rays at the geocenter. They differ from corresponding ICRS (a, 5) values because of 
annual aberration, annual parallax and gravitational light deflection due to the gravitational fields 
of the solar system bodies (apart from the Earth) and are independent of Earth's rotation. If these 
GCRS places, however, will ever play a role in practice is not clear. 

In the past apparent places of stars that were annually published e.g., in the 'Apparent Places 
of Fundamental Stars' played a role for certain problems. These places are related with the old 
traditional astronomical reference system, i.e., with some equator and equinox of date. Now with the 
ICRS we have a highly precise astronomical reference system that is basically independent of Earth's 
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rotation parameters and their determination. For several applications, however, the introduction 
of quantities like apparent places might still be useful especially if there is a reference to the local 
plumb line, i.e., to the zenith and the astronomical (or nautical) triangle can be employed. In that 
case Resolutions B1.7 and B1.8 (IAU 2001; Rickman 2001) come into play. These two Resolutions 
define some intermediate system that can be used for the definition of an intermediate position 
("inter, Winter) by the Celestial Intermediate Pole (CIP) and the Celestial Ephemeris Origin (CEO). 
Such intermediate position can be considered as modern version of the apparent place, defined in 
the GCRS. 

For astrometry at microarcsecond accuracies neither GCRS places nor intermediate places 
likely will play a role. To avoid problems related with non-linearities it is simpler to use an overall 
BCRS picture to describe not only the light-rays and the motion of gravitating bodies but also the 
trajectory of an observer. In that case only catalog and observed positions will be of importance. 



Resolution B1.3 concerns the definition of Barycentric Celestial Reference System (BCRS) and 
Geocentric Celestial Reference System (GCRS). The BCRS is defined with coordinates (ct, x l ) = 
where t = TCB. The BCRS is a particular version of the barycentric reference system of the solar 
system. The Resolution recommends to write the metric tensor of the BCRS in the form 



A comparison reveals that this form of the metric presents an extension of (2). Whereas the 
old form contains only, the Newtonian potential U the new one contains a scalar potential w and 
a vector potential w l . 

Actually the equations for goo and goi from (8) without the order symbols C(c -5 ) are always 



to the concrete form of the resolution we have added order symbols in (8). E.g., for goo the order 
symbol indicates that terms of order c -5 will systematically be neglected as stated in the notes to 
the Resolution. With these forms for goo an d goi one finds that spatially isotropic coordinates x l 
exist such that gij from equation (8) with the potential w from goo solves Einstein's field equations 
to first post-Newtonian order. Note that the form of (8) implies the barycentric spatial coordinates 
x l satisfy the harmonic gauge condition (see e.g. Brumberg & Kopeikin 1989; Damour, Soffel & 
Xu 1991). At this point, because of the freedom in the time coordinate, many different 'time gauge 



3.2. The Barycentric Celestial Reference System 




(8) 



correct and can be simply considered as definitions of w and w % in terms of goo and 90i- I n contrast 
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conditions' are still possible. The Resolution proceeds by recommending a specific kind of space 
and time harmonic gauge. One argument in favor of the harmonic gauge is that tremendous work 
on General Relativity has been done with the harmonic gauge that was found to be a useful and 
simplifying gauge for many kinds of applications. Moreover, the harmonic gauge condition (e.g., 
Weinberg 1972; Fock 1959) 



sT^ = 0, (9) 

where T^ u are the Christoffel-symbols of the metric tensor, is not restricted to some post-Newtonian 
approximation, but can be defined in Einstein's theory of gravity without any approximations. This 
may be important for future refinements of the IAU framework. With the harmonic gauge condition 
the post-Newtonian Einstein field equations take the form 



-\^ + A)w = -A7rGa + 0(c~ 4 ), (10) 



c 2 dt 2 

Aw i = -AirGa* + 0(c~ 2 ) . (11) 

Here a and a 1 are the gravitational mass and mass current density, respectively. Mathematically 
they are related to the energy-momentum tensor by 

a = -1 (T 00 + T ss ) , a* = -T \ (12) 

c z ' c 

The energy- momentum tensor T^ u generalizes the density p of the Poisson equation (3). In rela- 
tivity, energy density, pressure and stresses all act as source of the gravitational field. This implies 
that different kinds of energy contribute to the gravitational sources: kinetic energy, gravitational 
potential energy, energy of deformation etc., but since the kinetic energy depends upon the state of 
motion of the matter, the energy-momentum tensor that really acts as gravitational source shows 
a nontrivial transformation behaviour if we go from one reference system to another. In prac- 
tice, however, the energy-momentum tensor will usually not appear explicitly. This is because the 
gravitational potentials w and w l from (10)— (11) are completely determined by a and a 1 which 
can be considered as primary quantities. If we deal with problems where gravitational fields play 
a role only outside of astronomical bodies and admit a useful convergent expansion in terms of 
multipole moments (potential coefficients) only corresponding integral characteristics of the bodies 
like masses, quadrupole moments etc. show up explicitly, which are defined in terms of a and o % 
and whose numerical values will be fixed by observations. Because of (11) w % is sometimes called 
the gravitomagnetic potential since it results from mass currents (moving or rotating masses) just 
as the electromagnetic vector potential results from electric currents in Maxwell's theory of elec- 
tromagnetism. 
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Equation (10) generalizes the Poisson equation (3), hence the scalar potential w presents a 
relativistic generalization of the Newtonian potential U. Because of problems related with homoge- 
neous solutions and boundary conditions mathematically it is clear that these differential equations 
do not fix the harmonic solutions uniquely. Assuming space-time to be asymptotically flat (no 
gravitational fields far from the system), i.e., 

lim 9lxv = diag(-l, +1, +1, +1) (13) 

t— const 

the recommended solution reads 



«,(t,x) = G J M^^L + ^G^ J dVa(i,x')|x-x'|, (14) 

w%x) = G ! dV ^ (f,X 2 ■ (15) 
J |x — x'| 

It is obvious that the second time derivative term in (14) results from the corresponding operator in 
the field equation (10). This operator modifies the Laplacian from the Newtonian Poisson equation 
to the d'Alembertian and the similarity between the harmonic post-Newtonian field equations and 
Maxwell's equations of electromagnetism in the Lorentz gauge becomes obvious (actually one might 
replace the Laplacian by the d'Alembertian in (11) to post-Newtonian accuracy). From Maxwell's 
theory it is well known that the retarded potential solves the corresponding field equation 

w ret (t,x i ) = G [ rf 3 MWx') (16) 
J |x x I 

with / 

i ret = t - |X ~ X| . (17) 

C 

One might then expand the retarded potential in terms of 1/c. Note, that such an expansion also 
yields a term proportional to 1/c. If we stay within the first post-Newtonian approximation these 
1/c terms vanish due to the Newtonian mass conservation law. Such odd powers of 1/c indicate 
time asymmetric terms, i.e., they break time reversal symmetry. It is well known that such time 
asymmetric terms appear only to higher post-Newtonian order and will not be considered here. For 
that reason the retarded potential (16) leads to the recommended solution above. 

Comparing the form of the metric tensor in (8) with other forms that can be found in the 
literature (e.g., Will 1993) one might get the erroneous impression that something is missing in 
(8), which is not the case. If matter is discribed by some fluid model then formally (w,w l ) might 
be split into various pieces resulting from kinetic energy, gravitational potential energy, specific 
internal energy density, pressure etc. and the equivalence of our form of the metric tensor e.g. with 
that given in Will (1993) can be shown. This is explicitly demonstrated in Appendix B. 
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The point, however, is that a split of (a, a 1 ) of our metric potentials (w, w l ) or of the metric 
tensor itself into various pieces is usually unnecessary. If only gravitational fields outside the relevant 
bodies play a role (as is typically the case in celestial mechanics and astrometry) it is advantageous 
to keep such pieces together, since it will be the sum that determines the observables. One might 
argue, U is the 'Newtonian potential' and the rest can be identified as 'relativistic corrections'. 
This way of thinking, however, can be very misleading and presents a source of errors. As has 
been shown in the literature (e.g., Damour, Soffel & Xu 1991, 1993) suitably defined potential 
coefficients based upon w (not U) and w l can be introduced that can be determined from satellite 
data. From a more theoretical point of view the introduction of (w, w l ) has the advantage that the 
field equations (10)-(11) are formally linear, although the corresponding metric is not (because of 
the u> 2 -term). We used the word 'formally' since a depends upon w implicitly. This nonlinearity 
has been explicitly treated, e.g., by Brumberg & Kopeikin (1989), but this dependence becomes 
practically irrelevant if the fields outside of some matter distribution are parametrized by means of 
potential coefficients. This linearity implies that for an ensemble of iV-bodies 



N N 
w(t,x) = J2 w A(t,x), u/(i,x) = ^u4(t,x), (18) 
A=l A=l 

where the index A indicates the contribution related with body A where the integrals have to be 
taken over the support of body A only. This linearity, however, does not imply that body-body 
interaction terms have been neglected. If written explicitly wa will in general contain contributions 
from bodies B ^ A (see e.g., Eq. (54)). 

The BCRS metric tensor from the IAU 2000 Resolution B1.3 extends the form of the metric 
tensor given in the IAU 1991 Resolutions such that its accuracy is sufficient for most applications 
in the next years. Note that an extension of the old metric (2) is necessary (and has been in use 
for decades) for the derivation of the relativistic equations that form the basis of any modern solar 
system ephemeris (such as the JPL DE ephemerides) . Resolution B1.3 formalises this extension. 



3.3. The Geocentic Celestial Reference System 

Resolution B1.3 continues to define the GCRS which represents a particular version of the local 
geocentric reference system for the Earth. Its spatial coordinates X a are kinematically non-rotating 
with respect to the barycentric ones (see e.g., Brumberg & Kopeikin 1989; Klioner & Soffel 1998). 
The geocentric coordinates are denoted by (T, X), where T = TCG. In the relation between x % 
and X a from Resolution B1.3 let us replace the unit matrix 5 a i by a general rotation matrix R a i 



X a = R n 



1 



+ 0(c- 4 ), 



where = x — x#. If the two sets of spatial coordinates are aligned for all times, i.e., if R a 



as 
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is the case for the GCRS spatial coordinates, then X a is defined to be kinematically non-rotating 
with respect to the barycentric spatial coordinates x l . The Resolution recommends writing the 
metric tensor of the GCRS in the same form as the barycentric one but with potentials W(T, X) 
and W a (T,X). Explicitly 

1 2W 2W 2 ml _ 5n 
G 00 = -l + ^ ^ + °( c )- 



G 0a = -^W a + O(c- 5 ), (19) 
G ab = S ab (l + %w) +0(c- 4 ) 



and the geocentric field equations formally look the same as the barycentric ones (10)-(11) but 
with all variables referred to the GCRS. Again one decisive advantage of this recommendation is 
the formal linearity of the field equations. This linearity admits a unique split of the geocentric 
metric into a part coming from the Earth itself and a remaining part resulting from inertial and 
tidal forces. Therefore it is recommended to split the potentials W and W a according to 



W{T, X) = W E (T, X) + W cxt (T, X), W a (T, X) = W E (T, X) + W^(T, X). (20) 

The Earth's potentials We and W E are defined in the same way as we and w E , (i.e., equations 
(14)-(15) with integrals taken over the volume of the whole Earth) but with quantities calculated 
in the GCRS. Outside the Earth the potentials (W, W a ) admit a power series expansion in terms 
of R = |X| and all negative powers of R are contained in We and W E . For that reason the Earth's 
potentials admit multipole expansions that look very similar to the Newtonian ones. This point 
will be discussed below in more detail. 

It is useful to split the external potentials W cxt and W^ xt further. They can be written in the 

form 



W cxt = W" tida i + W-ner, = W^al + ^iner , (21) 

where the tidal terms are at least quadratic in X a and the inertial contributions Winer and W^ er 
are just linear in X a . Explicitly, 



Winer — Qa X a , 

W-^ er = -\c 2 e abc 9!i ncr X c . (22) 

Mathematically the Q a term is related with the 4-acceleration of the geocenter in the external 
gravitational field, a quantity that vanishes for a purely spherical and non-rotating Earth (for a 
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mass monopole more precisely) that moves along a geodesic in the external gravitational field. The 
Q a term therefore results from the coupling of higher order multipole moments of the Earth to the 
external tidal gravitational fields (to the external curvature tensor of space-time in mathematical 
language) . Q a characterizes the deviation of the actual worldline of the origin of the GCRS from 
a geodesic in the external gravitational field. With 

Wext(*,x) = ^ w A(t,*) , «4rt(*> x ) = ^2 ™a(*> x ) 

to Newtonian order Q a is given by 

Qa = 5ai (Jl~[ w cxt (x E ) - a l ^j . (23) 

Here, x l E (t),v E (t) = dx E /dt and a E = dv E /dt are the barycentric coordinate position, velocity 
and acceleration of the origin of the GCRS (geocenter). The appearance of 5 a i results from the 
fact that the GCRS is defined as kinematically non-rotating with respect to the BCRS. The reason 
to retain 5 a i in the transformations here and below results from the desire to distinguish between 
BCRS (spatial indices taken from the second part of the latin alphabet, starting with the letter i) 
and GCRS quantities (spatial indices taken from the first part of the alphabet). 

The full post-Newtonian expression for Q a (denoted by G a (T) in the Damour-Soffel-Xu papers) 
can be derived from (6.30a) of Damour, Soffel & Xu (1991). To get an idea about orders of 
magnitude the absolute value of Q a due to the action of the Moon is of order 4 x lCr n m/s 2 
(Kopeikin 1991). 

The term W^ er describes a relativistic Coriolis force due to the rotation of the GCRS with 
respect to a dynamically non-rotating geocentric reference system. Such a rotation has several 
components, often referred to as geodesic, Lense-Thirring and Thomas precessions 

fWr = flop + ^ltp + ^tp (24) 

with 



3 _ 

2 

^ltp = — ^Vxw cxt (x B ), (25) 
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in obvious notation. As a relativistic precession, the geodesic precession, flap, is proportional 
to 1/c 2 . It is also proportional to the bary centric coordinate velocity ve and the gradient of 
the external gravitational scalar potential u> ex t at the geocenter (the barycentric coordinate ac- 
celeration of the geocenter to sufficient accuracy). The order of magnitude is given by |ficp|~ 
1.5(ve/c)(GM s /c 2 AU)(c/AU) ~ 3 x 10~ 15 s _1 ~ 2"/century. Thomas precession is also propor- 
tional to 1/c 2 and the barycentric coordinate velocity of the geocenter but also to the geodesic 
deviation term Q a . The order of magnitude of Thomas precession is |fixp| ~ 0.5(ve/c)\Q\/c ~ 
7 x 10~ 24 s _1 ~ 4 x 10~ 9 "/century, i.e., negligible with respect to geodesic precession. 

Finally the Lense-Thirring precession results from the gradient of the external gravitomagnetic 
potential at the geocenter. If we consider some spherically symmetric solar system body A, then 
the gravitomagnetic potential W A of it is given by (see (49) below) 

G {8a x X)" 
Wa ~ 2 IF 

in its own local rest frame. Transformation into the BCRS according to the rule indicated below 
in (31) leads to 



w A (t,x)=G 



(S A xr A y | M A i 
2r5 + r A VA 



where = x - and is the barycentric velocity of body A. In our case the spin and motion 
of our Sun and Moon will give the dominant contributions to fiurp: |^ltp| ~2x lCT 3// /century. 

The definition of the GCRS implies that the spatial GCRS coordinates X are kinematically 
non-rotating with respect to the BCRS ones, x (as indicated by the (5 a ;-term in Resolution B1.3). 
Because of geodesic precession locally inertial coordinates precess with respect to the GCRS by an 
amount of | Oj ner | = 1.9198"/century (Brumberg et al. 1992). Let us forget about the mass of the 
Earth and imagine a torque free gyroscope at the geocenter, moving along the actual trajectory 
of the geocenter. It will precess by this amount in our GCRS. Since the GCRS does not present 
a locally inertial reference system Coriolis forces caused by geodesic precession-nutation appear 
in every GCRS dynamical equation of motion, e.g., of Earth's satellites. As recommended by 
IERS Conventions (2003) these additional forces should be taken into account. Moreover, geodesic 
precession-nutation has to be considered in the precession-nutation model formulated in the GCRS. 
E.g., the basic post-Newtonian equation of Earth's intrinsic angular momentum S reads 



— + ft incr X S = D , (26) 

where D is the external torque (Damour, Soffel & Xu 1993). As long as observations of Earth's 
orientation parameters are referred to the GCRS they will contain geodesic precession-nutation 
automatically. 
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Because of the eccentricity of the Earth orbit the leading term in ficp has an annual and a 
semi-annual part that leads to geodesic nutation in longitude with 

AV> G p = 0. 153 sin I' + 0.002 sin 21' , (27) 

where the amplitudes are in mas and I' is the mean anomaly of the Earth-Moon barycenter 
(Fukushima 1991; Brumberg et al. 1992; Bois & Vokrouhlicky 1995). 

Wtidai * s a g eneranza ti° n of the Newtonian tidal potential 

W^T n (T, X) = w cxt (x E + X) - w cxt (x E ) - X • V Wcxt (x E ) . (28) 

Full post-Newtonian expressions for Wtidai and VF/L j can be found in Damour, Soffel & Xu (1992). 
There lf ex t is denoted by W and a tidal expansion in powers of local spatial coordinates by means 
of suitably defined tidal moments is given in (4.15) of Damour, Soffel & Xu (1992). Expressions 
for Wtidai and W t " dal in closed form are given in (Klioner k, Voinov 1993). The quadratic term as 
dominant term of Wtidai reads 



Wtidal 



1=2 



^G^ al X a X b . (29) 



If the external bodies are taken as mass monopoles the explicit expression for G*^ al (not to be 
confused with the GCRS metric tensor) is given in (3.23) of Damour, Soffel & Xu (1994). Higher- 
order terms in this approximation can be found in Klioner, Soffel, Xu, Wu (2000). 

Finally, the local gravitational potentials We and W E of the Earth are related to the barycen- 
tric gravitational potentials we and w l E by (5f = 5 l a = 8 a j) 

W E (T 1 X)=w E (t,x) + -|^^(t,x) + 0(c- 4 ), 

W E (T,X) = 5f («4(t,x) -v E w E (t,x)) +0(c- 2 ) (30) 

or by the inverse transformation 



w E {t,-x.) = W E (T,X) (l + ^vl) +^S ta v E W E (T,X) + 0(c- 4 ), 

w E (t, X ) = 5l W E (T, X) + v% W E (T, X) + 0(c~ 2 ). (31) 

The relations between the geocentric gravitational potentials W and W a , and the barycentric ones 
w and w l follow from the coordinate transformations between the BCRS and GCRS discussed 
below. 
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3.4. Coordinate transformations 

The metric tensors in the BCRS and GCRS allow one to derive the rules for the transformations 
between the BCRS coordinates x M and the GCRS ones X a from the tensorial transformation rules. 
It is obvious that these transformations can be written in two equivalent forms: i) as x M (T, X a ) or ii) 
as X a (t, x l ). Whereas the first form was used in the Damour-Soffel-Xu formalism (Damour, Soffel 
& Xu 1991, 1992, 1993, 1994), the second one was presented in the Brumberg-Kopeikin formalism 
(Brumberg & Kopeikin 1989; Kopeikin 1988; Brumberg 1991; Klioner & Voinov 1993). It should 
be pointed out that the transformation from one version to the other is not so trivial because of the 
barycentric coordinate position of the geocenter that appears in the first form as function of TCG 
and as function of TCB in the second one. In Resolution B1.3 T = TCG and X a are presented 
as functions of t = TCB and x l . The explicit form of the transformations is given in the text 
of Resolution B1.3 (see, Appendix A below). Apart from the terms of order (9(|X| 3 ) that appear 
in the time transformation of order 0(c -4 ) all terms can be obtained from the results derived by 
Kopeikin (1988) and Damour, Soffel & Xu (1991). The cubic and higher order terms in |X| as 
represented by the function C have been derived by Kopeikin (1988) and analyzed in full detail by 
Klioner & Voinov (1993). As is also clear from Klioner & Voinov (1993) the expression for C is 
not unique, but only constrained by the gauge and field equations so that the simplest possibility is 
an expression for C containing cubic terms only. It is this simplest expression that is recommended 
in Resolution B1.3. 

The full 4-dimensional coordinate transformation is just an extension of the usual Lorentz 
transformation. Indeed, if we neglect all gravitational fields and acceleration terms then the coordi- 
nate transformation in Resolution B1.3 can be written in the form (r = x — xe(£), (3 = v/c = const.) 



T=t ( 1 - \f ~ \f) ~ (l + \?) ^ + 0(c-% 
X = r+I(v.r)J + 0( C - 4 ). (32) 



If we now write X£(i) = vi we obtain 



X = x- (l + i/3 2 )v( + i(vx) J +0(c- 1 ). (33) 

which is nothing but a Lorentz tranformation from Special Relativity Theory 

(v • x\ (W — 1) 

*--2-J ' X = X -^+V (V ' X)V 
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in the corresponding approximation since 

7 ee (1 - /? 2 )-V2 = 1 + \f + 1(3* + 0(c" 6 ) . 

Note that the inverse transformations are obtained simply by replacing (t, x) by (T, X) and the 
velocity v by —v. 

Neglecting the 1/c 4 terms in the T — t relation given in Resolution B1.3 one gets 

T = t -^([ t { ( ! Y + ™ext(x E )) dt + v%r%j + 0(c- A ) (34) 

which reduces to the old recommendation (4) since t = TCB, T = TCG and w cxt (^-E(t)) reduces 
to U ex t(t, xb(£)) in the Newtonian limit. The more accurate version of this transformation will be 
discussed below. 

Let us also note that the BCRS, GCRS and the transformation between them have been 
discussed by Klioner & Soffel (2000) in the framework of the PPN formalism with parameters (3 
and 7. For the limit of General Relativity (3 = 7 = 1 all the formulas given in that publication 
become equal to those derived in the framework of the new IAU Resolutions that refer solely to 
Einstein's theory of gravity. 

3.5. Potential coefficients 

3.5.1. General post- Newtonian multipole moments 

For many problems it is advantageous to present the local gravitational potentials of the 
Earth as multipole series that usually converge everywhere outside the Earth. To this end one 
has to introduce a certain set of multipole moments or potential coefficients for the Earth. A 
certain set of potential coefficients, called Blanchet-Damour (BD) moments (Blanchet & Damour 
1986; Damour, Soffel Sz Xu 1991) defined to first post-Newtonian order has especially attractive 
features. Moreover, by using such Blanchet-Damour moments we get a very simple form of the 
multipole expansion of the post-Newtonian potentials (these expansions have almost Newtonian 
form). Basically two sets of BD moments occur in the formalism: mass multipole moments and 
spin multipole moments. Theoretically these moments can be derived from the distribution of mass 
and matter currents inside the body but for an observer they simply present parameters which can 
be directly estimated from observations. 

Expressed in terms of symmetric and trace-free Cartesian tensors the BD moments are denoted 
by Ml and Sl- Here, L is a multi-index of I different indices all taking the values 1,2,3, i.e., 
L = i\%2 ...ii and every index i = 1, 2, 3. Explicit expressions for Ml and Sl as integrals over the 
Earth can be found in, e.g., Blanchet & Damour (1986); Damour, Soffel & Xu (1991) 
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M L {T)= / d 6 XX^ + 



d 2 



2(2/ + 3)c 2 dT 2 
4(2/ + 1) d 



(/ + l)(2/ + 3)c 2 cff 



/ 



d 3 XX L X 2 T l 
d 3 XX aL Y, a 



I > 0, 



(35) 



Sl(T)i 



d^Xe 



ab<ci ^"L— l>ayjfe 



/ > 1, 



(36) 



where the integrations extend over the body under consideration and 



e(t,x) = - (r 00 + r ss ), 



£ a (T,X) = -T 



Oa 



(37) 



Here T Miy = T^ U {T, X a ) are the components of the energy-momentum tensor in the GCRS. Both 
the caret and the sharp brackets indicate the symmetric and trace-free (STF) part of the object 
or of the indices enclosed by the brackets (see, e.g. Damour, Soffel & Xu (1991) p. 3277 for the 
explicit definition of the STF part of an object). Some basic information on the operations with 
the STF objects can be found, e.g. in Blanchet & Damour (1986); Damour, Soffel & Xu (1992). 

For practical applications, however, their explicit form will not be needed, since these quantities 
are parameters characterizing the gravitational field of the corresponding body which are fitted to 
observations. The set Ml is equivalent to a set of potential coefficients C[ m and Si m that appear in 
a much more familiar spherical harmonic expansion of We- The first non- vanishing spin moment 
(the spin dipole) of a body agrees with its spin vector (total intrinsic angular momentum). The 
multipole expansion of We and W% reads 



M L d L ^- + ^M L d L \x\ 



+ 4tAt + 0(c- 4 ), 



(38) 



oo . 



\i r 



MaL-idh-iT^r + j^-^e abc S C L-id bL _i—— 



K a + 0{c- 2 ), (39) 



where 



oo 



^(Z + l)! 2Z + 3 



2l + 1 v » 1 
1X1 ' 



(40) 
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V L = [ S a X aL d 3 X. (41) 
Jv 

Here the dot stands for d/dT and 8l for d l jdx %1 . . . dx 1 '. Also the comma denotes partial differen- 
tiation, A >T = and A >a = dA/dX a . 

The gauge function A does not enter the post-Newtonian equations of motion. The latter 
contains only the BD multipole moments Ml and Sl- The only place where the function A should 
be accounted for is in the transformation between the various time scales. However, these gauge 
terms are of order 1/c 4 in the metric tensor so for the problem of clock rates they are basically 
of second post-Newtonian order. These terms are much less than 10~ 18 in the geocentric metric 
tensor and will be neglected. For that reason the A-terms are not mentioned in Resolution B1.4. 

3.5.2. Approximate expansion of the scalar gravitational potential 
A spherical harmonic expansion of We equivalent to (38) without the A term reads (R = |X|) 



W E (T,X) 



GM 



E 



R 

+0(c- 4 ) 



oo +Z 
1=2 m=0 



Re 
~R 



Pl m (cos 9) (Ci m (T, R) cos Tficj) -\- 5; m (T, R) sin mc, 



(42) 



with 



CF m {T,R) = Ct{T) - ^L_|-J^C&(T) (43) 

St(T,R) = St{T) - 2(^1)5^2^^ • ^ 

Let us stress that as stated in Resolution B1.4 C^(T) and Sf m (T) refer to the GCRS coordinates 
and are related with approximately constant potential coemcents in a terrestrial system that is 
rotating with the Earth (i.e. those from an Earth's model) by time-dependent transformations. 
For a rigid axially symmetric body rotating about its symmetry axis with angular velocity the 
second time derivative terms will vanish. Let us estimate these terms for the Earth. From the 
order of magnitude of the I = m = 2 terms in the reference system corotating with the Earth one 
finds C§2 an d 5^ of order 1CT 6 . The expected order of magnitude of the second time derivative 
terms is (fig Re/c) 2 ~ 10~ 12 times smaller than the corresponding 'Newtonian terms' from or 
S§2- The Newtonian terms lead to contributions in Gqq of order 10 -15 and hence the second time 
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derivative terms to contributions of order 10 . This is about nine orders of magnitude less than 
the 2W 2 /c 4 term in Goo which is of order 10~ 18 . For that reason these second time derivative terms 
in the Earth's metric can safely be neglected at present. They are not mentioned in Resolution 
B1.4. 



3. 5. 3. Approximate expansion of the vector gravitational potential 

Let us now come to the gravitomagnetic vector potential of the Earth, W^. As can be seen 
from (39) this potential is determined by the set of spin moments and the first time derivative of the 
mass moments. As already mentioned, to characterize the gravitational field outside of some matter 
distribution in GRT two independent sets of multipole moments have to be used that in principle 
should be determined from observational data. So far the spin moments of some astronomical body 
have not been studied and more work is needed here. Formally the spin moments of the Earth are 
given by expression (36) above. Since for the post-Newtonian metric we need these spin moments 
only to Newtonian order we might proceed with a simple Newtonian model of a rigidly rotating 
Earth with 

s = e (n x x) , 

where £ is the gravitational mass-energy density in the GCRS and f2 is the angular velocity of 
rotation that at this place has to be defined only to Newtonian order. Under this assumption all 
spin moments are proportional to the angular velocity and one might define a set of Cartesian 
tensors such that 

S L = C Ld n d . (45) 
These tensors obey the following Newtonian relations 

C Ld = -M Ld + ill^JV^ , (46) 

where 

M L = f Sl L d 3 X , N L = [ SX 2 1 L d 3 X . (47) 
Je Je 

Note that CLd is symmetric and trace-free only in the first L indices. Moreover, for the Newtonian 
mass moments Ml one has 



M L = ~C <L> ■ 



(48) 
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For a homogeneous (E = const.) and spherical Earth with radius Re one finds for I = 1 the usual 
expression for the moment of inertia tensor 

Cab = Sab {^ MR l) 

that yields the total intrinsic angular momentum (spin) vector of the Earth according to S a = 
C a b Q h ■ For a spherically symmetric and mass centered Earth all mass moments Ml with / > 1 
vanish and also all quantities Nl with I > 0. Hence in such a simple model only the spin vector 
is different from zero and all higher spin moments vanish. For that reason we also considered a 
rigidly rotating homogeneous oblate spheroid with equatorial radius A and polar radius C. For 
such a model all even spin moments vanish since they are proportional to Cl with odd I. On the 
contrary, odd spin moments proportional to Cl with even I are nonzero. For the spin dipole the 
usual result Cxx = Cyy = M(A 2 + C 2 )/5 and Czz = 2MA 2 /5 for the moment of inertia tensor 
is found. By means of computer algebra all components Cl can be found for any value of I. Let 
rj = (4MA 4 /525)e 2 with e 2 = (A 2 - C 2 )/A 2 ~ 2f, where / is the usual flattening. Assuming 
f2 d = (0, 0, fi) we found all non-vanishing I = 3 terms, up to symmetries and terms of order f 2 : 
Sxxz = Syyz = 3r]Q and Szzz = — 6??f2. This implies that the metric term resulting from the 
spin octupole of the Earth near the Earth's surface is about 10 4 times smaller than the one from 
the spin dipole. In the following the contributions of higher spin moments will be neglected. 

Besides the spin moments the first time derivative of the mass moments contribute to the 
gravitomagnetic field of the Earth. For I = we encounter a A4 a term that vanishes if the post- 
Newtonian center of mass condition M a = is imposed. The next term is given by M a b that 
is of order \Cf 2 \ MR 2 E £l and would vanish for an axially symmetric rigid body rotating about its 
symmetry axis as well as the time derivative of all higher mass moments. For the Earth the 
M. a b term is smaller than the spin term (which is of order 2MR 2 E Q,/b) by a factor determined 
by C§2 — 1-6 x 10~ 6 and hence negligible. On the other hand, the vector potential W E (T, X) is 
employed only in the calculations of small relativistic effects (e.g. Lense-Thirring effects, higher- 
order relativistic effects in the time transformations). This implies that the expansion (39) for 
W E (T, X) can be truncated to the approximate expression 

^(r,x) = -f£^, (49) 

where is a vector with components <S a . This expression can be found in many standard textbooks 
on GRT (Weinberg 1972; Will 1993) and is usually related with Lense-Thirring effects resulting 
from the Earth's rotational motion. 

The reason why to characterize W E by the spin vector and not by the angular velocity vector 
of the Earth is a conceptual one since usually it is advantageous to characterize the gravitational 
field of the Earth in the outside region by multipole moments. To get W E the Earth's spin vector 



-27- 



is needed only to Newtonian order and can be taken from current precession-nutation models. 
Although one might use Newtonian concepts to relate the gravitomagnetic field of the Earth with 
some Earth's angular velocity, we prefer to employ the well-defined concept of multipole moments 
here which are independent of any theoretical assumptions on the rotational motion of the Earth. 



In the gravitational N-body problem the potential coefficents of a body A are defined in its 
corresponding local reference system (analogous to the GCRS for the Earth) . For many applications 
it is sufficient to keep only the mass monopoles of the solar system bodies, i.e., to put 



for all bodies and to keep the masses only, i.e., each body A is characterized by the value for its 
post-Newtonian mass A4a (we also put Pl = 0). In the following we will use the symbol Ma 
instead of Ma to be compatible with the text of the IAU 2000 Resolutions. 

From the transformation rules for the metric potentials (31), expansions (38)-(39) and formula 
(18) one derives the metric in the barycentric coordinate system in the form (8) with 



3.6. The barycentric metric in the mass monopole approximation 



M L = for 1>1, 



S L = for I > 1 



(50) 



w = w Q 5-A 



(51) 



where 



w (t,x) = 



GM A 



(52) 



A 



and 



A(*,x) = £A A (t,x) 



(53) 



A 



with {yba = Xfj — x.A and a^ = d\A/dt) 




(54) 
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Furthermore, in our approximation 

W %x)=Y;^v a . (55) 

Note, that we have chosen the minus sign in front of A to have a plus sign in the c~ 4 part of goo (see 
Resolution B1.5). Note furthermore, that the post-Newtonian Einstein-Infeld-Hoffmann equations 
of motion for a system of mass monopoles, that form the basis of modern solar system ephemerides, 
can be derived from that form of the barycentric metric (for details see Damour, Soffel & Xu 1991). 
Thus, the barycentric mass monopole metric given above is already in use for the description of 
the solar system dynamics. 

One improvement of this simple mass monopole model is to consider the spin dipoles of the 
various bodies as well (that is, to consider also to be non-zero). Actually Resolution B1.5 is 
based upon such a mass monopole spin dipole model where modifications from the simple mass 
monopole model are indicated explicitly. 



4. Time and frequency applications in the solar system 

For practical applications concerning time and frequency measurements in the solar system 
it is necessary to consider a conventional model for the realization of time coordinates and time 
transformations. This model should be chosen so that i) its accuracy is significantly better than the 
expected performance of clocks and time transfer techniques, ii) it is consistent with the general 
framework of Section 3 and iii) it may readily be used with existing astrometric quantities, e.g. 
solar systems ephemerides. 

Regarding item (i), we may derive reasonable accuracy limits for such a model in a straight- 
forward way. At present the best accuracies are reached by Cs-fountain clocks operating at less 
than two parts in 10 15 in fractional frequency (Lemonde et al. 2001; Weyers et al. 2001). Their 
frequency stability for time spans up to a few days characterized by a standard Allan deviation 
is of order a y { T ) = 4 x 10- 14 r~ 1 / 2 , for an integration time r in seconds. In the near future, high 
accuracy laser cooled Rb clocks (Bize et al. 1999) and space-borne Cs clocks (Lemonde et al. 
2001) are expected to reach accuracies of a few parts in 10 17 in fractional frequency and stabilities 
of order a y (r) = 1 x 10~ 14 r -1 / 2 . The uncertainty in the time transformations should induce errors 
that are always lower than the expected performance of these future clocks. Including a factor 
2 as safety margin, we therefore conclude that time coordinates and time transformations should 
be realized with an uncertainty not larger than 5 x 10~ 18 in rate or, for quasi-periodic terms, not 
larger than 5 x 10~ 18 in rate and 0.2 ps in amplitude. 

For the spatial domain of validity of the transformations, we note that projects like SORT 
(Solar Orbit Relativity Test) plan to fly highly accurate clocks to within 0.25 AU of the Sun, 
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which is therefore the lower limit for the distance to the barycenter that we will consider. In the 
geocentric system we will consider locations from the Earth's surface up to geostationary orbits 
(|X| < 50 000 km). 

To comply with item (ii) , we render the developments following the general framework outlined 
in Section 3 and we show (hi) how the time transformations, e.g., TCB — TCG, may be performed 
with the existing astrometric quantities and tools. 

4.1. Barycentric reference system 

Let us write the barycentric metric potential w(t, x) in the form 

w = w + w L - 4tA , (56) 

where wl contains contributions from higher potential coefficients with I > 1 and can be determined 
from equation (38) and the transformation rules of the metric potentials. Evaluating the A^ terms 
from Resolution B1.5 (equation (54) plus spin terms) for all bodies of the solar system, we find 
that in the metric tensor |A^(f,x)|/c may reach at most a few parts in 10 17 in the vicinity of 
Jupiter and about 10~ 17 close to the Earth. Presently, however, for all planets except the Earth, 
the magnitude of A^(i,x)/c 4 in the vicinity of the planet is smaller than the uncertainty in wq/c 2 
or wl/c 2 originating from the uncertainties in its mass multipole moments so that it is practically 
not needed to account for these terms. Nevertheless, when new astrometric observations allow 
to derive the moments with adequate uncertainty, it will be necessary to do so. In any case, for 
the vicinity of a given body A, only the effect of A^(t,x) is needed in practice, i.e., the effect of 
X^b^a ^s(i,x) is smaller than our accuracy specifications. For the comparison of proper time of 
a clock in the vicinity of the Earth with that of other clocks in the solar system or with TCB it 
may thus be needed to account for A#(£,x)/c 4 . 

From (8) and (56) the transformation between proper time of some observer and TCB may 
be derived within our accuracy limit 



dr/dTCB = 1 - 1 (w + w L + + 1 V - ^v 2 w + W + ±v$ + a) , (57) 

where v l is the BCRS coordinate velocity of the observer. Similarly, the transformation between 
TCB and TCG in the immediate vicinity of the Earth, accurate to the limits as specified above can 
be derived from the general post-Newtonian TCB — TCG transformation from Resolution B1.3 as 



TCB -TCG = c 



s 



+ u>o,cxt(x£) dt + v l E r l E 
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3 1 

^v 2 E w , cxt (^E) + 4v l E w l ext (x E ) + 2^o,cxt( x s) ) dt 



J t 

-(3u;o,ext(x E ) + v 2 E /2)v E r E ], (58) 



where t is TCB. Here u>o,ext is defined by (52) with summation over all solar system bodies expect 
the Earth. Note that to was not explicitly defined in Resolution B1.5 (2000). It is the origin of 
TCB and TCG, defined in Resolution A4 (1991) (see Section 2.1). The external contributions to 
wl and A are beyond our accuracy limit and can be neglected here. 

This equation is composed of terms evaluated at the geocenter (the two integrals) and of 
position dependent terms linear in \r E \, terms with higher powers of \r E \ having been found to be 
negligible. The integrals may be computed from existing planetary ephemerides (Fukushima 1995; 
Irwin & Fukushima 1999). Since, in general, the planetary ephemerides are expressed in terms of 
a time argument T cph = (1-L B ) TCB + T c ° ph (Standish 1998; Irwin & Fukushima 1999), the first 
integral will be computed as 



+ w ,ext(x£;) dt 



- cph 



epho 



+ W0,ext(x£) cff cph 



/(I -L B ). 



(59) 



Terms in the second integral of (58) are secular and quasi periodic. They amount to ~ 
1.1 x 10~ 16 in rate {dTC B / dTCG) and primarily a yearly term of ~ 30 ps in amplitude (i.e. 
corresponding to periodic rate variations of amplitude ~ 6 x 10~ 18 ). Position-dependent terms 
in c~ 4 are not negligible and reach, for example, an amplitude of 0.4 ps (~ 3 x 10~ 17 in rate) in 
geostationary orbit. 



4.2. Geocentric reference system 

Evaluating the contributions of the different terms in the metric tensor of the GCRS given in 
Resolution B1.3 to the dr/d(TCG) transformation on the Earth's surface and up to geostationary 
orbit we find that terms of order c~ 2 reach 7 parts in 10 10 while the contributions from W 2 and W a 
do not exceed 5 parts in 10 19 . Also, the terms from W- mcv in W remain below 2 x 10 -20 . Therefore, 
the terms given in the IAU 1991 framework with the metric of the form (2) are sufficient for time 
and frequency applications in the GCRS in the region |Xj < 50 000 km for present and foreseeable 
future clock accuracies. Note, some care needs to be taken when evaluating the potential W at 
the location of the clock which is not trivial when accuracies of order 10 -18 are required (Klioner 
1992; Petit & Wolf 1994; Wolf & Petit 1995). 

Presently, the time scale of reference for all practical issues on Earth is Terrestrial Time (TT) 
or one of the scales realizing it and differing by some time offset (e.g., TAI, UTC, GPS time). TT 
was defined in IAU Resolution A4 (1991) as: "a time scale differing from the Geocentric Coordinate 
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Time TCG by a constant rate, the unit of measurement of TT being chosen so that it agrees with 
the SI second on the geoid" . According to the transformation between proper and coordinate time, 
this constant rate is given by d(TT) /d(TCG) = 1 — U g /c 2 = 1 — Lc, where U g is the gravity 
(gravitational + rotational) potential on the geoid (this notation is used instead of the usual Wo' 
to avoid confusion with GCRS gravitational potential W used throughout the paper). 

Some shortcomings appeared in this definition of TT when considering accuracies below 10~ 17 . 
First, the uncertainty in the determination of U g is of order 1 m 2 s -2 or slightly better (Bursa 1995; 
Groten 1999). Second, even if it is expected that the uncertainty in U g improves with time the 
surface of the geoid is difficult to realize (so that it is difficult to determine the potential difference 
between the geoid and the location of a clock). Third, the geoid is, in principle, variable with 
time. Therefore it was decided to desociate the definition of TT from the geoid while maintaining 
continuity with the previous definition. The constant Lq was turned into a defining constant 
with its value fixed to 6.969290134 x 10~ 10 (Resolution B1.9, see Appendix A). This removes the 
limitations mentioned above when realizing TT from clocks on board of terrestrial satellites (such 
as in the GPS). In Table 1 we present numerical values for the constants Lc, Lc and Lb relating 
the mean rates of different relativistic timescales. 



5. Final remarks 

The relativistic Resolutions of the IAU represent a post-Newtonian framework allowing one to 
model any kind of astronomical observations in a rigorous, self-consistent manner with accuracies 
that are sufficient for the next decades. They replace the old IAU relativistic framework that was 
insufficient for many reasons discussed above. These new Resolutions, however, are not expected 
to lead to dramatic changes. In fact in many fields of application the models presently in use 
are already fully compatible with the new IAU Resolutions and in this sense the IAU Resolutions 
officially fix the status quo. Let us give some examples for that. 

The metric tensor of the BCRS allows one to derive the Einstein-Infeld-Hoffman equations 

IAU 1991 IAU 2000 IAU 2000 

s/s s/s ms/yr 



L c 1.480813 x 10~ 8 1.48082686741 x 10~ 8 467.313 

L G 6.969291 x 10~ 10 6.969290134 x 10~ 10 21.993 

L B = L C + L G - L c L G 1.550505 x 10~ 8 1.55051976772 x 10" 8 489.307 



Table 1: The constants relating the mean rates of different relativistic timescales. Both the val- 
ues adopted by the IAU 1991 recommendations and the IAU 2000 resolutions are given. As an 
illustration the IAU 2000 values are also given in milliseconds per Julian year. 
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which have been used since the 70s to construct the JPL numerical ephemerides of planetary motion 
(Newhall, Standish & Williams 1983). The BCRS is the basic astrometric reference system where 
concepts such as 'radial velocity' or 'proper motion' are defined (Lindegren & Dravins 2003). The 
metric tensors of both GCRS and BCRS and the transformation between corresponding coordinates 
were used to formulate the VLBI model employed by the IERS since 1992 (IERS Standards 1992; 
IERS Conventions 1996, 2003). The equations of motion of the Earth's satellites recommended 
by the IERS (IERS Standards 1992; IERS Conventions 1996, 2003) are compatible with the new 
IAU framework and can be derived from the given metric tensor of the GCRS. 

The models used for costructing the HIPPARCOS catalog make it clear that this catalog 
represents a materialization of the BCRS. The full power of the new IAU theoretical framework 
will be needed to construct the model for astrometric positional observations with an accuracy of 
1 microarcsecond which will be necessary for future astrometric missions. Work in this direction 
has already started (Klioner 2003). 

It is obvious that this explanatory supplement presents only a first step to show how the new 
IAU Resolutions concerning relativity should be employed in practice. Much more work will be 
necessary to reach that goal. 

The anonymous referee is thanked for his valuable suggestions to improve the text and to make 
it more readable. 
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A. IAU Resolutions concerning Relativity Adopted at the 24th General Assembly 
Resolution B1.3: Definition of Barycentric Celestial Reference System and 



The XXVIth International Astronomical Union General Assembly 
Considering 



1. that the Resolution A4 of the XXIst General Assembly (1991) has defined a system of space- 
time coordinates for (a) the solar system (now called the Barycentric Celestial Reference 
System, (BCRS)) and (b) the Earth (now called the Geocentric Celestial Reference System 
(GCRS)), within the framework of General Relativity, 

2. the desire to write the metric tensors both in the BCRS and in the GCRS in a compact and 
self-consistent form, and 

3. the fact that considerable work in General Relativity has been done using the harmonic gauge 
that was found to be a useful and simplifying gauge for many kinds of applications, 

Recommends 

1 . the choice of harmonic coordinates both for the barycentric and for the geocentric reference 
systems, 

2. writing the time-time component and the space-space component of the barycentric metric 
g^u with barycentric coordinates (t, x) (t = Barycentric Coordinate Time (TCB)) with a 
single scalar potential w(t,x) that generalizes the Newtonian potential, and the space-time 
component with a vector potential w l (t, x); as a boundary condition it is assumed that these 
two potentials vanish far from the solar system, 



Geocentric Celestial Reference System 



explicitly, 



2w 2w 2 
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with 



w(t,*) = G J d 3 x'^l + ^G^ J dV<x(t,x')|x-x'| 
(t,x) = G J , 



„,<(/. x) - r; / d 3 x'^^- 



Here, <r and are the gravitational mass and current densities, respectively. 

3. writing the geocentric metric tensor G a p with geocentric coordinates (T, X) (T = Geocentric 
Coordinate Time (TCG)) in the same form as the barycentric one but with potentials W(T, X) 
and W a (T, X); these geocentric potentials should be split into two parts — potentials We 
and W E arising from the gravitational action of the Earth and external parts W ext and W® xt 
due to tidal and inertial effects; the external parts of the metric potentials are assumed to 
vanish at the geocenter and admit an expansion into positive powers of X, 

explicitly, 



2W 2W 2 
G 00 = -l + -^ 

G 0a = ~^W a , 



c 3 



The potentials W and W a should be split according to 



W(T, X) = W E (T, X) + W cxt (T, X), 
W a (T, X) = W E (T, X) + W? xt (T, X). 



The Earth's potentials We and W E are defined in the same way as w and w a but with 
quantities calculated in the GCRS with integrals taken over the whole Earth. 

4. using, if accuracy requires, the full post-Newtonian coordinate transformation between the 
BCRS and the GCRS as induced by the form of the corresponding metric tensors, 

explicitly, for the kinematically non-rotating GCRS (T = TCG, t = TCB, r E = x i - x E (t) 
and a summation from 1 to 3 over equal indices is implied), 
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T = t 



X a = 6 ai 



A(t) + 44 
i fi 



+ 



Bit) + B\ty E + + c(i, x) 



+ 0(c" 5 ), 



c 2 V2 



+ 0(c- 4 ), 



where 



A ( t ) = n V E + ^cxt(xE), 



13 1 

^B(t) = --v% - -v 2 E w cxt (x E ) + 44^(xe) + -^ xt ( XE ), 

B i {t) = --v 2 E v E + 4«4 t (x E ) - 3^w C xt(x E ), 

d d 1 

B ij (t) = - V y aj Q a + 2^— tf)j xt (x E ) - 4^jW«t(3CE) + 2^cxt(x E ), 

C(t,x) = -lr|(d^). 

Here x^, v^, and are the barycentric position, velocity and acceleration vectors of the 
Earth, the dot stands for the total derivative with respect to t, and 



d 



^-Wext(x E ) - a l E 



Q a = S ( 

The external potentials, u> ex t and w l ext , are given by 



A+E A^E 



where E stands for the Earth and wa and w\ are determined by the expressions for w and 
w l with integrals taken over body A only. 



Notes 



It is to be understood that these expressions for w and w % give goo correct up to 0(c~ 5 ), goi up to 
0(c~ 5 ), and gij up to 0(c~ A ). The densities a and a 1 are determined by the components of the 
energy momentum tensor of the matter composing the solar system bodies as given in the references. 
Accuracies for G a p in terms of c~ n correspond to those of g^ v . 
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The external potentials W cxt and W® xt can be written in the form 

W ex t=W / tidal + W iner , 

Wtidai generalizes the Newtonian expression for the tidal potential. Post- Newtonian expressions for 
Wtidai an d ^tidai can ^ e found in the references. The potentials Wi ncT , W^ er are inertial contri- 
butions that are linear in X a . The former is determined mainly by the coupling of the Earth's 
nonsphericity to the external potential. In the kinematically non-rotating Geocentric Celestial Ref- 
erence System, Wf nci describes the Coriolis force induced mainly by geodesic precession. 

Finally, the local gravitational potentials W E and W E of the Earth are related to the barycentric 
gravitational potentials we and w E by 

W E (T,X)=w E (t,x) (l + ^v E ^j -^v E w E (t,x) + 0(c~ 4 ), 
W a E (T, X) = 6 ai (w E (t, x) - v E w E (t, x)) + 0(c- 2 ) . 
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Resolution B1.4: Post-Newtonian Potential Coefficients 

The XXVIth International Astronomical Union General Assembly, 

Considering 
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1. that for many applications in the fields of celestial mechanics and astrometry a suitable 
parametrization of the metric potentials (or multipole moments) outside the massive solar- 
system bodies in the form of expansions in terms of potential coefficients are extremely useful, 
and 

2. that physically meaningful post-Newtonian potential coefficients can be derived from the 
literature, 



Recommends 



1. expansion of the post-Newtonian potential of the Earth in the Geocentric Celestial Reference 
System (GCRS) outside the Earth in the form 



W E (T,X) 



GM e 
R 



1 + E E (if J P ^ cos 9 )( C L(T) cosm0 + Sf m (T) sinmt 



1=2 m=0 



Here and are, to sufficient accuracy, equivalent to the post-Newtonian multipole 
moments introduced by Damour et al. (Damour et al., Phys. Rev. D, 43, 3273, 1991). 9 and 
4> are the polar angles corresponding to the spatial coordinates X a of the GCRS and R = \X\, 
and 

2. expression of the vector potential outside the Earth, leading to the well-known Lense-Thirring 
effect, in terms of the Earth's total angular momentum vector Se in the form 



W E (T,X) 



G (X x S E ) a 



Resolution B1.5: Extended relativistic framework for time transformations and 
realization of coordinate times in the solar system 



The XXVIth International Astronomical Union General Assembly, 



Considering 



1. that the Resolution A4 of the XXIst General Assembly (1991) has defined systems of space- 
time coordinates for the solar system (Barycentric Reference System) and for the Earth 
(Geocentric Reference System), within the framework of General Relativity, 
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2. that Resolution B1.3 entitled "Definition of Bary centric Celestial Reference System and Geo- 
centric Celestial Reference System" has renamed these systems the Barycentric Celestial 
Reference System (BCRS) and the Geocentric Celestial Reference System (GCRS), respec- 
tively, and has specified a general framework for expressing their metric tensor and defining 
coordinate transformations at the first post-Newtonian level, 

3. that, based on the anticipated performance of atomic clocks, future time and frequency mea- 
surements will require practical application of this framework in the BCRS, and 

4. that theoretical work requiring such expansions has already been performed, 
Recommends 



that for applications that concern time transformations and realization of coordinate times within 
the solar system, Resolution B1.3 be applied as follows: 

1. the metric tensor be expressed as 



goo z 
9iy 



1 - -o(«>o(t,x) + w L (t,x)) + -r(u;§(t,x) + A(t,x)) 



f «;*(*, x), 



1 + 



2w (t,x) 



where (t = Barycentric Coordinate Time (TCB), x) are the barycentric coordinates, wq = 
G J2a Ma/ta with the summation carried out over all solar system bodies A, ta = x-x^, x^ 
are the coordinates of the center of mass of body A, ta = |i"a|, and where wl contains 
the expansion in terms of multipole moments [see their definition in the Resolution B1.4 
entitled "Post-Newtonian Potential Coefficients"] required for each body. The vector potential 
w l (t,-x) = ^^(^x) an d the function A(i,x) = J2a Aa(*,x) are given in note 2. 

2. the relation between TCB and Geocentric Coordinate Time (TCG) can be expressed to 
sufficient accuracy by 



TCB - TCG = c 



S t (if + Wo ' ext ( X£ )) dt + 

j t (~8~^ ~ ^^0,ext(x B ) + 4v E W l cxt (x E ) + ^0,ext( x s)^i 



3w , 



cxt 



>e) + ^-)v e t e 



where v E is the barycentric velocity of the Earth and where the index ext refers to summation 
over all bodies except the Earth. 
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Notes 



1. This formulation will provide an uncertainty not larger than 5 x 10 in rate and, for quasi- 
periodic terms, not larger than 5 x 10™ 18 in rate amplitude and 0.2 ps in phase amplitude, for 
locations farther than a few solar radii from the Sun. The same uncertainty also applies to the 
transformation between TCB and TCG for locations within 50000 km of the Earth. Uncertainties 
in the values of astronomical quantities may induce larger errors in the formulas. 

2. Within the above mentioned uncertainties, it is sufficient to express the vector potential w l A (t,x) 
of body A as 

-(r A xS A y | M A v A ~ 
2r\ r A y 

where S A is the total angular momentum of body A and v A is the barycentric coordinate velocity of 
body A. As for the function A A (t, x) it is sufficient to express it as 



w A (t,x) = G 



A A (i,x) = 



GM A 

r A 



-K + £ 

B+A 



GMb 
TBA 



k„,k\2 



+ 



(r\v 



A U A) 



+ r k A a k A 



+ 



2Gv k A (r A x S A ) 



where tba = |x_b — and a A is the barycentric coordinate acceleration of body A. In these 
formulas, the terms in S A are needed only for Jupiter (S ~ 6.9 x 10 38 m?s~ 1 kg) and Saturn (S ~ 
1.4 x 10 38 m 2 s~ 1 kg), in the immediate vicinity of these planets. 



3. Because the present recommendation provides an extension of the IAU 1991 recommendations 
valid at the full first post-Newtonian level, the constants Lq and L b that were introduced in the 
IAU 1991 recommendations should be defined as < TCG /TCB >= 1 - L c and < TT/TCB >= 
1 — Lb, where TT refers to Terrestrial Time and <> refers to a sufficiently long average taken 
at the geocenter. The most recent estimate of Lc is (Irwin, A. and Fukushima, T., 1999, Astron. 
Astroph., 348, 642-652) 

L c = 1.48082686741 x 10~ 8 ±2x 10~ 17 . 

From Resolution B1.9 on "Redefinition of Terrestrial Time TT", one infers Lb = 1.55051976772 
xl0~ 8 ±2 x 10 -17 by using the relation 1 — Lb = (1 — Lq)(1 — Lc). Lc is defined in Resolution 
B1.9. 
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Because no unambiguous definition may be provided for Lb and Lq, these constants should 
not be used in formulating time transformations when it would require knowing their value with an 
uncertainty of order 1 x 10~ 16 or less. 

4- IfTCB — TCG is computed using planetary ephemerides which are expressed in terms of a time 
argument (noted T cp ^) which is close to Barycentric Dynamical Time (TDB), rather than in terms 
ofTCB, the first integral in Recommendation 2 above may be computed as 



The XXVIth International Astronomical Union General Assembly, 
Considering 



1. that IAU Resolution A4 (1991) has defined Terrestrial Time (TT) in its Recommendation 4, 



2. that the intricacy and temporal changes inherent to the definition and realization of the geoid 
are a source of uncertainty in the definition and realization of TT, which may become, in the 
near future, the dominant source of uncertainty in realizing TT from atomic clocks, 

Recommends 



that TT be a time scale differing from TCG by a constant rate: dTT/dTCG = 1 — Lq, where 
L G = 6.969290134 x 1(T 10 is a defining constant, 

Note 




-f+w , ex t(xE)jdt /(l -L B ). 



Resolution B1.9: Re-definition of Terrestrial Time TT 



and 



Lq was defined by the IAU Resolution A4 (1991) in its Recommendation 4 as equal to Uq/c 2 where 
Uq is the geopotential at the geoid. Lq is now used as a defining constant. 
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B. Comparison of the IAU metric with versions given in the literature 

In this Appendix we will compare the metric (8) with well-known results from the literature. 
Only for this purpose we will consider the material composing the various bodies of the system to 
behave like an ideal fluid (for the IAU 2000 Resolutions this is not assumed). In the ideal fluid case 
the energy-momentum tensor can be written in the form 



Here, p denotes the rest-mass density, p is the pressure, II is the specific internal energy (e.g., Will 
1993), v l (t, x) is the velocity of the corresponding material element and 




r 0i = pci; i + O(c- 1 ), 

T ij = pv i v j + P 5 ij + 0{c- 2 ). 



(Bl) 




(B2) 



From (12) and (Bl) we derive 




a^pv 1 + 0(c~ 2 ). 



(B3) 



Introducing the metric potentials 




(B4) 



(B5) 



we obtain from (14) and (15) 



w = U + 2$i + 2<J> 2 + $3 + 3<I>4 



w l = V t + 0(c~ 2 ), 



(B6) 



with 
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Vi = j ^ dV , X = ~ G j P'\* ~ x 'l ^ ' ( B7 ) 
and the comma denotes partial differentiation as in 

The post-Newtonian metric (in harmonic gauge) can then be written as 

500 = -1 + ^ [2U + 4*i + 4$ 2 + 2$ 3 + 6$ 4 ] - ^U 2 + , 

90i = ~Vi, (B8) 

gij = Sij I 1 + ^2" 

To compare, e.g., with the metric in (Will 1993) we transform from harmonic coordinates, used in 
the present paper and recommended by the IAU, to standard post-Newtonian (SPN) coordinates 
used by several authors including Will. This is achieved by a gauge transformation of the form 
(e.g., Eq. (3.12) of Damour, Soffel & Xu (1991)), 

wspn = w - 4j A, t , ^spn = w l + (B9) 

with 

X = -\x,t- (BIO) 

This implies that the x-term disappears from w and hence from <7 o when standard post-Newtonian 
coordinates are employed, but the goi term gets affected by this transformation. Using the relation 



X,« = Vi-Wi 

with 



^'.(x-x-)]^-**)^ 



one verifies that the metric induced by the potentials (14) and (15) agrees in General Relativity 
with (5.28) of (Will 1993). 
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